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We report on our simulation of finite temperature lattice QCD with two flavors of O(a) Symanzik-improved 
fermions and 0(a 2 ) Symanzik-improved glue. Our thermodynamic simulations were performed on an 8 3 x 4 
lattice, and we have performed complementary zero temperature simulations on an 8 3 x 16 lattice. We compare 
our results to those from simulations with two flavors of Wilson fermions and discuss the improvement resulting 
from use of the improved action. 
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1. INTRODUCTION 

The study of finite temperature QCD with 
Wilson-type quarks is desirable in order to esti- 
mate any systematic errors of similar simulations 
with Kogut-Susskind quarks. However, Wilson 
thermodynamics has proved to be difficult and 
burdened with lattice artifacts Q|. It is plausible 
that an action which converges to the continuum 
action faster in the a — > limit would be cured 
of such spurious effects. 

2. ACTION 

For the gauge action, we start with the one 
loop, on-shcll Symanzik improved action derived 
by Liischer and Weisz ||. We implement the 
tadpole improvement scheme in order that lat- 
tice perturbation theory be more convergent |||| . 
We choose to define the "mean link" uq and the 



strong coupling constant a s through the plaque- 
tte §-§: 

The coefficients of the rectangle operator and 
the twisted 6-link operator, /3 rec t and Awist re- 
spectively, are given in terms of the coefficient of 
the plaquette (3 and uq as in M: 
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In practice, we estimate uq in a self-consistent 
manner: we tune it so that it agrees with the 
fourth root of the space-like plaquettes. 

The Wilson fermion action has errors of 0(a). 
The Symanzik improvement program is used to 
improve the action ||. After tadpole improve- 
ment the fermion action is 
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where Sw is the usual Wilson fermion action, and 
iFp V is the familiar clover-shaped link operator. 
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Figure 2. Polyakov loop vs. hopping parameter 
for 8 3 x 4 improved Wilson thermodynamics. 



Figure 1. Phase diagram of Symanzik-improved 
action. Octagons represent the Nt = 4 ther- 
mal crossover, and diamonds indicate estimates 
of vanishing pion mass. Zero temperature simu- 
lations were performed at the crosses. 



3. RESULTS 

Our thermodynamics simulations were done on 
an 8 3 x 4 lattice at six fixed values of (3 while 
varying k across the thermal crossover (tuning 
uq self-consistently at each parameter set). We 
used the hybrid Monte Carlo algorithm and col- 
lected data from at least 1000 trajectories for the 
simulations in the crossover region. Furthermore, 
zero temperature simulations on an 8 3 x 16 lattice 
were performed in order to provide hadron masses 
in the region of the thermal crossover line. The 
phase diagram (figure |l|) summarizes our run pa- 
rameters. 

Figure || shows the Polyakov loop as a function 
of k for the six values of (3. One can observe 
that the transition appears steeper for stronger 
coupling: a feature also present in Nt — 4 Wilson 
thermodynamics R|. Still, the crossover for the 
improved action does not appear to be as steep 
as for the unimproved action. 



One would like to make direct comparison be- 
tween the two actions of their respective crossover 
behavior without depending on the bare param- 
eters. In this work, we use measurements of the 
lattice pion mass squared at values of k near the 
crossover. Then, we can plausibly overlay curves 
of thermodynamic observables for two actions run 
at comparable m n /mp. Below we list m^/mp 
along the Nt = 4 crossover for both clover and 
Wilson [pi actions. 
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Using measurements of the pion mass near the 
crossover region Q , we can interpolate in order to 
estimate («m,) 2 as a function of 1/k. Then, we 
can plot the thermodynamic observables against 
the pion mass squared. This shows that the 
crossover is indeed smoother for Nt = 4 clover 
than N t = 4 Wilson (see figure |J) . 

Finally, the confinement-deconfinement tem- 
peratures for different two-flavor lattice actions 
are shown in figure ^|. One consequence of our im- 
provement scheme is to lower the Wilson Nt = 4 
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Figure 3. Polyakov loop vs. pion mass squared. 
The crossovers (N t = 4) for both actions, at the 
couplings shown, occur at the same n — p mass 
ratio: m n /m p = 0.83. 



critical temperature at a given mass ratio. This 
brings the calculation of T c /m p into better agree- 
ment with Nt = 6 Wilson and with staggered 
fermion thermodynamics. There is one important 
caveat: we do not yet know if the clover simula- 
tions have reached a plateau in /m p . If T c /m p 
continues to rise at lower m^/rrip (lower f3) then 
the aforementioned agreement is accidental. Fur- 
thermore, we should remember that T c /m p tends 
toward zero as m w /m p approaches unity, i.e. as 
m q — >■ oo. Measurements of the string tension will 
provide a scale which is insensitive to the quark 
mass. 

4. CONCLUSIONS 

We have shown that the Nt = 4 thermal 
crossover is smoother for the Symanzik-improved 
action. It could be that T c /m p is in better 
agreement with staggered fermion results; how- 
ever running at the thermal crossover at lower 
m^/rrip is needed to confirm this. 

This work was supported by the U.S. Depart- 
ment of Energy and the National Science Foun- 



Figure 4. Crossover temperature in units of the 
p mass vs. m^/rrip. 

dation. 
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